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Introduction 


1 . 1  Summary  1 

The  Department  of  Meteorology  and  Wind  Energy  at  the  Riso 
National  Laboratory  in  Denmark  participated  and  supported: 


The  MADONA: 

(Meteorological  And  Diffusion  Over  Non-uniform  Area) 

Trials 


Conducted  at 
Porton  Down,  Salisbury 
CBDE-UK;  Sep.  7-25  1992. 


Riso  contributed  with  the  following  activities: 


1)  Continuous  smoke  dissemination  from  ground  level  sources. 

2)  High-resolution  protile-measurements  of  instantaneous 
plume  dispersion  and  concentration  fluctuations  by  use  of 
a  mini-LIDAR  system.  (Remote  sensing) 

3)  Micro-meteorological  measurements  of  wind  and  temperature 
quantities  from  two  7  meter  tall  met-towers  equipped  with 
20  Hz  sonic  anemometers  for  characterization  of  the  mean 
and  turbulent  state  of  the  boundary  layer. 

4)  Mean  wind  and  diffusion  modelling  using  our  real-time 
diagnostic  models  LINCOM/RIMPUFF . 


1.2  Objectives 


The  scientific  objective  for  our  participation  in  MADONA  was  to 
investigate  fundamental  questions  regarding  the  statistical 
properties  of  the  instantaneous  concentrat ion  fluctuations  in 
smoke  plumes  and  puffs.  The  joint-venture  Porton  Down  ex¬ 
periments,  with  international  participation  and  support 
facilitated  this  objective.  Modelling  objectives  included 
investigation  and  experimental  evaluation  of  flow  and  diffusion 
modelling  based  on  our  LINCOM/RIMPUFF  real-time  codes. 

With  our  high-resolution  mini-lidar  system,  sequential  measure¬ 
ments  (1/3  Hz)  were  obtained  of  the  instantaneous  cross  wind 
concentration  profiles  at  fixed  down  wind  distances  from  the 
source.  The  source  consisted  of  ground  level  continuous  smoke 
releases  and  elevated  smoke  puff  releases. 

The  cross  wind  concentration  profiles  were  measured  with  an 
effective  spatial  resolution  of  1.5  meter. 

From  the  conducted  puff  and  plume  diffusion  experiments,  we  have 
by  post-processing  and  data  reduction  obtained  important  con¬ 
centration  fluctuation  statistics  including: 

Mean-  and  variance  profiles,  intermittency  profiles, 
probability  distribution  functions,  dispersion  parameters 
(including  instantaneous  puff  width).  In  addition  we  dis¬ 
tinguish  statistics  obtained  in  a  fixed  frame  (absolute 
diffusion) ,  and  center-of-mass  frame  (relative  diffusion) . 


Derived  statistics  (to  be  performed)  include: 

Dispersion  coefficients,  and  their  dependence  on  sampling 
time  (from  instantaneous  to  length  of  experiments,  typi¬ 
cally  1-hour. ) 

Two-particle  distance-neighbor  functions  yielding  infor¬ 
mation  about  the  instantaneous  spatial  correlation  func¬ 
tion  of  in-plume  fluctuations. 

Extreme  statistics  of  exeedance  of  certain  reference 
levels,  and  the  (ensemble  mean)  duration  of  such  exeedan- 
ces . 


Our  Lidar  based  concentration  fluctuation  measurements  were 
supported  by  our  micro-meteorological  observations  from  two 
tower  mounted  sonic  anemometers,  in  addition  to  the  joint  mean 
wind  and  turbulence  data  gathered  simultaneously  by  the  partici¬ 
pating  groups  at  Porton  Down. 


1.3  Background: 


Within  the  CEC  Radiation  Protection  Research  Programme  for 
Nuclear  Safety,  Riso  is  commissioned  to  conduct  reference  smoke 
plume  dispersion  experiments,  where  we  based  on  the  mini-LIDAR 
facility  are  studying  concentartion  fluctuations  for  dispersion 
model  evaluation  and  model  uncertainty  assessments.  Also,  we  are 
conducting  full  scare  reference  and  training  experimental  data 
for  the  CEC  supported  real-time  decision  support  system  RODOS, 
where  atmospheric  dispersion  is  a  key  element. 

The  layout  and  experimental  set-up  for  the  Porton  Down  MADONA 
experiments,  including  the  many  other  participants,  provided  an 
unique  opportunity  to  combine  and  jointly  perform  already 
scheduled  uncertainty  and  evaluation  experiments  with  the 
general  objectives  for  MADONA. 

Our  entire  scientific  outcome  from  the  smoke  diffusion  trial  are 
consequently  offered  for  inclusion  in  the  MADONA  data  base, 
including  the  basic  scientific  analysis  as  described  above. 

All  raw  and  statistical  analyzed  data  from  our  LIDAR  and  micro 
meteorological  measurements  are  shared.  The  latter  includes  10- 
min  mean  and  turbulent  scaling  parameters  such  as  variances, 
energies,  shear  stress  and  heat  flux  (Monin-Obukhov  length). 

1.4  Reporting: 

This  is  a  first  interim  data  report  describing  the  experimental 
setups,  time  periods  and  atmospheric  conditions  for  the  MADONA 
trials . 

A  final  report,  including  the  statistical  analysis  and  inter¬ 
pretation  of  results  will  be  available  end  of  august  1993. 

1.5  Acknowledgements: 

Riso  obtained  additional  financial  support  for  expenses 
associated  with  moving  our  experimental  facilities  and  personnel 
to  Porton  Down, UK. 

Our  participation  in  the  MADONA  trials  were  hereafter  made 
possible  by  support  from  the  CEC  Radiation  Protection  Research 
Programme,  DGXII,  under  contract  No.  BI7-0017-C,  and  by  a  grant 
from  the  US  Army  Research  Branch,  London,  under  contract 
No .  DAJA45-92-M-0344  .  Mr.  Ron  Cionco  at  the  U.S.  Army  Atmospheric 
Science  Laboratory,  White  Sands,  MN. ,  USA  and  Mr.  Jerry  Comati 
at  the  US  Army  Research  Branch,  London,  are  gracefully  ack¬ 
nowledged  . 
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Technical  describtion 


2 ■ 1  Mean  and  turbulence  measurements  from  two  micro-meteoro¬ 
logical  masts: 

Two  7  meter  high  meteorology  mast  were  erected  and  instrumented 
at  the  Porton  Down  site  near  the  Westerly  source  point,  see  Fig. 
II. 1.  (Porton  Down  landscape  seen  from  South-West.) 

The  net-towers  were  each  equipped  with  Kaijo  Denki  DAT  300  sonic 
anemometers  in  order  to  obtain  both  mean-profile  and  turbulence 
scaling  parameters  for  the  MADONA  dispersion  trials. 

Appendix  II. 1:  Contains  directly  measured  mean  and  turbulence 
quantities,  including:!.)  Mean  Wind  speed, 2) 
Direction, 3)  Turbulent  Kinetic  Energy,  4) 
Stress,  5)  Sensible  Heatflux,  6)  Monin-Obukhov 
Stability  parameters  and  7)  Temperature  statis¬ 
tics  . 

Appendix  II. 2  Comparison  between  5  and  10  min  averaged  statis¬ 
tics  for  day:  23/9  192;  sonic  1  and  sonic  2. 

The  measurements  were  obtained  from  two  acoustic  three-dimen¬ 
sional  sonic-anemometer/thermometers  (Kaijo  Denki  type  DAT -300) , 
located  at  7  meter  level  above  the  floor. 

Running  continuously  during  the  entire  experimental  campaign, 
they  provided  us  with  real-time,  on-line  surface  layer  tur¬ 
bulence  statistics  of  the  most  important  scaling  parameters  for 
dispersion.  By  use  of  our  PC-based  on-line  data  acquisition 
system  and  associated  computers,  we  measured  and  calculated  (in 
real-time)  the  entire  co-variance  matrix  of  the  fluctuating  qu¬ 
antities  u'jV’jW1  and  T'  (by  the  correlation-method),  and  from 
here,  the  velocity  and  temperature  intensities,  the  shear-stress 
and  the  sensible  heat  flux  was  deduced  at  10-min  intervals  in 
real  time.  With  the  additional  information  of  the  inversion  hei¬ 
ght  from  rawin-sonde  balloons  (not  provided  by  Rise) ,  the  sonic 
anemometer  data  enables  us  to  calculate  the  very  important 
scaling  parameter  for  convective  conditions:  W  .. 

Our  sonic-based  on-line  "turbulence  monitoring-system"  provides 
the  experimenters  and  the  project-manager  with  real-time 
measurements  of  the  atmospheric  stability  and  turbulence  state 
at  Porton  Down,  concurrently  with  each  diffusion  test  and 
provided  also  direct  measurements  of  the  variances  and  fluxes 
requested  for  the  subsequent  model  evaluation  study. 

The  P.iso  Sonic-DAS-system  was  put  in  continuous  operation 
throughout  the  whole  experimental  period. Our  directly  measured, 
10-min  averaged  quantities  included : <u 1 u •> ,  <u'v'>,  <u'w'>, 

<u'T’>,  <v'v'>,  <v'w’>,  <v'T'>,  <w ’ w ' > , <w ' T ' > ,  <T'T’>,  in 

addition  to  the  mean  wind  speed  and  mean  temperature. 

The  accuracy  of  the  sonic  anemometer/thermometer  type  DAT  300 
is,  according  to  the  manufacturer  (Kaijo  Denki  Co.):  +/-  l%  on 
measured  velocity  and  temperature  signals,  the  resolution  is  0.5 
cm/s  and  0.025  K,  respectively.  It  measures  all  signals  20  times 
per  second  and  has  a  10  Hz  low-pass  output  filter.  For  mea¬ 
surements  taken  7  meter  above  the  ground,  this  instrument  will 


provide  very  adequate  resolution  of  the  turbulent  signals,  both 
with  respect  to  spatial  and  temporal  resolution:  The  on-line 
data  acquisition  system  (DAC_SYS)  point-samples  the  four  signals 
U,V,W  and  T  at  a  rate  of  10  Hz,  and  subsequently  calculates  the 
co-variance  matrix  after  each  10-min  period  based  on  the 
corresponding  6000  readings  of  each  signal. 

The  directly  measured  time  series  (consisting  of  four  times  6000 
measurements  each  per  10-min  period)  are  also  stored  for 
subsequent  spectral  analysis. 

Data-,  both  mean  and  of  the  turbulence  quantities,  have  already 
been  processed,  stored  and  backed-up  as  10-min  average  quanti¬ 
ties.  Additional  to  visual  inspection  of  raw  time  series,  more 
ob}ective  quality  assurance  in  accordance  with  Hojstrup,  1993 
will  be  performed  on  central  data  periods. 


2 . 2  High-resolution  Time-series  and  Spectral  Analysis: 

Additional  spectral  analysis  can  at  request  be  performed  on  10- 
min  average  meteorological  quantities.  For  each  of  the  scheduled 
diffusion  tests,  the  following  is  available  on  request: 

1.  Time  series  plots  of  high-resolution  (10  Hz)  wind  and  tem¬ 
perature  signals  from  the  sonic  anemometer. 

2.  Time  series  of  short-time  averaged  (1-min  running  mean) 
quantities  of  wind  speed,  direction,  turbulence  level, 
shear  stress  and  heat  flux. 

3.  Spectral  analysis  of  the  three  wind  (u',v',w')  and  tempe¬ 
rature  (T1)  signals. 

This  information  is  of  importance  in  a  subsequent  modelling  eva¬ 
luation.  in  particular  for: 

I  Gaining  insight  into  the  fine  scale  temporal  evaluation  of 
the  boundary  flow  and  turbulence  in  the  valley  during  each 
diffusion  test. 

II  gaining  insight  into  which  turbulent  scales  are  of  impor¬ 
tance  and  most  energetic,  and  consequently  responsible  for 
the  exchange  and  diffusion  processes,  and 

III  Providing  detailed  turbulence  measurements  at  the  source 
point  for  use  in  the  subsequently  planned  model  simulation 
of  the  experiments. 

We  provide  on  request,  for  each  of  the  successfully  performed 
diffusion  tests: 

re  i,  10-Hz  digitized  raw  data  of  U,V,W  and  T  on  PC-dis- 

kette . 

re  ii.  Time  series  plots  of  sonic  data  and  of  the 
correspond-  ding  1-min  running  mean  turbulence 
quantities . 

re  iii.  Plots  of  the  three  velocity  spectra  and  the  tempera¬ 
ture  spectrum  f^om  0.001  to  5  H7 
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LIDAR-raeasurements  of  surface-released  smoke  plumes: 


The  Rise  National  Laboratory  operates  a  mini-LIDAR  system  for 
determination  of  instantaneous  concentration  profiles  inside 
dispersing  aerosol-plumes. 

The  system  yields  high  spatial  and  temporal  resolution  measure¬ 
ments  of  in-plume  concentration  fluctuations,  from  which 
important  and  basic  scientific  plume  dispersion  characteristics 
can  be  inferred. 

For  the  MADONA  trials,  we  provided  a  powerful  and  sturdy  ar¬ 
tificial  smoke  generator,  which  produced  ground-level  releases 
of  smoke  plumes  or  elongated  puffs. 

The  smoke  consisted  of  non-toxic  white  plume  of  sub-micron 
aerosols.  The  aerosols  were  conglomerated  SiO  2  and  NH  ^OH  that 
is  detectable  by  our  1.064  p  lidar  system.  In  this  way, 
measurements  of  the  cross-wind  near-ground  concentration 
profiles  were  obtained  at  various  distances  downwind. 

During  evenings  and  nights  with  stable  atmospheric  conditions, 
the  lidar  system  can  pick  up  the  smoke  plume  concentrations  out 
to  several  kilometers  from  the  source  point. 

At  Porton  Down,  the  mini-LIDAR  system  measured  instantaneous 
"snapshots"  every  3.  or  5.  second  of  the  aerosol  back-scatter 
cross-section . 

Table  III.l:  contains  information  of  lidar  positions  and 

measurement  periods  obtained  at  the  MADONA 
trials . 

Includes  also  drawings  of  source  and  lidar 
positions  for  the  total  of  23  experiments. 

Fig.  III.l:  LIu/rR  mean  concentrations  from  MADONA 

experiment  No.  22f  (1  hour  22  min  continuous 

release) . 


Fig.  III.2a-d:  From  Madona  exo.  No.  15h:  Lidar  measured  mean, 
variance,  fluctuation  intensity  and  intermitten- 
cy  crosswind  profiles  in  fixed  frame  (meander 
included) . 

Fig.  III.3a-d:  From  Madona  exp.  No.  15h:  Lidar  measured  mean, 
variance,  fluctuation  intensity  and  intermitten- 
cy  crosswind  profiles  in  moving  frame  (meander 
removed) . 


The  measurements  related  to  the  number  of  aerosols  in  the  small 
measurement  volume,  which  is  a  "ciqar-shaped"  cylinder  of  ap¬ 
prox.  1  meter  in  diameter  and  approx.  1.5  meter  long. 

Extensive  data-processing  of  the  raw  1 idar-signals  is  in 
progress  in  order  to  produce  accurate  and  quality-assured  cross- 
wind  aerosol  concentration  profiles.  In  return,  we  obtain  from 
each  experiment,  not  only  the  mean  concentration  profile,  but 
also  a  measure  of  the  natural  fluctuation  and  inherent  uncer¬ 
tainty  associated  with  dispersion.  This  is  of  importance  for 
subsequent  model  evaluation  studies. 

During  madona,  we  obtained  both  stable,  neutral  and  unstable 
experiments  based  on  both  plume  and  puff  releases.  The  release 
rate  of  the  aerosol  generator  was  kept  constant  during  in¬ 
dividual  experiments. 

Our  contribution  to  the  MADONA  data  base  contains,  for  each 
experiment : 

1)  Time  series  of  cross-wind  profiles 

2)  Plots  of  mean  cross-wind  concentration 

3)  Plots  of  fluctuating  (rms) 

4)  Plots  of  in-plume  intermittency 

5)  Measurements  of  the  entire  measured  concentration  pdf 
(probability  density  distribution)  from  which,  for  instan¬ 
ce,  different  exceedance  statistics  and  high-level  expo¬ 
sure  events  can  be  inferred. 

The  posterior  data  processing,  reduction  and  analysis  of  our 
Lidar-act ivity  is  rather  involved,  and  will  require  between  6 
ana  12  month  of  data  preparation  and  analysis  following  the 
experiments  to  finalize. 


Modelling 


4 • 1  Mean-flow  ncdelli ** a  directions : 


In  connection  with  our  on-going  activities  with  the  CEC-RODOS 
activities  for  real-time  on-line  decision  support,  we  ad¬ 
ditionally  setup  our  modelling  system  LINCOM/RIMPUFF  at  the 
MADONA  si*-': 

For  mean-flow  wind  field  calculations  at  the  Porton  Down 
experimental  area  ,  we  have  setup  our  fast  PC-version  of  LINCOM 
to  calculate  winds  over  the  experiment  site  at  various  heights 
and  for  different  synoptic  wind  directions.  LINCOM  is  a  non-hyd¬ 
rostatic  spectral,  neutral  wind  code,  cf.  Santabarbara  (1993). 


LINCOM  is  setup  to  run  based  on  the  measured  tower  winds  and  as 
such,  it  is  suitable  as  driver  for  dispersion  studies  over 
complex  terrain.  Alternatively,  the  model  can  run  on  a  single 
wind  speed  and  direction  as  input,  if  the  dominating  winds  are 
driven  by  synoptic  forces. 

Our  studies  include  comparison  of  measured  tower  winds  with 
model  results  based  on  single  inputs,  the  MADONA  study  with 
several  wind  observations  are  unique  for  this  model/data 
intercomparison . 

At  present  this  flow  model  is  under  development  in  order  to 
incorporate  the  effect  of  thermally  driven  forcing  such  as 
night-time  drainage  flows  and  day-time  up-valley  breezes.  Since 
this  model  is  a  central  as  a  "driver"  for  our  puf f-dispersion 
model  RIMPUFF  in  complex  terrain,  it  is  worthwhile  to  test  its 
applicability  for  the  Porton  Down  terrain  in  question,-  and 
possibly'  even  use  the  obtained  results  from  these  experiments 
in  orcier  to  evaluate  our  efforts  regarding  improving  the  model 
to  include  the  effects  of  thermal  forcing. 


4 . 2  Combined  Mean-wind  and  Dispersion  Models  for  Real-time 
Plume  Dispersion  Assessment  over  undulating  Terrain: 


Our  dispersion  modelling  work  is  centered  about  a  Lagrangian 
Mesoscale  Puff-diffusion  model  (RIMPUFF) ,  which  over  the  years 
has  obtained  acknowledgement  and  use  internationally,  for 
instance  in  the  UFOMOD  and  the  COSYME  consequence  assessment 
systems  available  in  public  domain  by  the  CEC.  For  dispersion 
calculations  over  complex  terrain,  the  PC-based  RIMPUFF  model 
is  presently  twinned  with  the  non-hydrostatic  mean  flow  model 
LINCOM,  which  is  based  on  the  Navier-Stokes  equations  for 
motion,  including  mass,  momentum  and  heat  conservation.  The 
combination  of  the  flow  and  diffusion  model  results  in  an 
extremely  fast  computer  code,  which  advantageously  can  be  tested 
and  evaluated  in  real  time  with  experimental  data  from  the 
present  diffusion  experiments. 

A  modelling  effort  based  on  selected  MADONA  diffusion  trials  is 


consequently 

in  progress. 

Appendix  IV: 

Contains  the  first  preliminary  model  results  of 
wind  and  diffusion  calculations  from  the  MADONA 
trials:  Sep  17, release  time  11:10  -  16:10. 

Fig. 

IV.  1: 

Shows  LINCOM  output  of  hourly  consecutive  wind 
field  calculations  based  on  tower  data  inputs. 

Fig. 

IV.  2  : 

Shows  detailed  stream  line  plot  of  one  of  these 
wind  fields,  including  tower  wind  data??? 

Fig. 

IV.  3  : 

Shows  instantaneous  plume  positions  at  ten 
minute  intervals  and  also  total  integrated 
concentration  (dose) . 

Fig. 

IV.  4  : 

Porton  Down  terrain  vertical  magnified  x  20 
seen  fj.om  south-west 

Fig. 

IV.  5: 

Example  of  tower  data  representation  for  Lincom 
inputs . 

Fig . 

IV.  6 : 

Example  of  Lincom  grid  wind  field  (10  meters 
height)  300  meter  horizontal  resolution 

Fig. 

IV.  7  : 

As  Fig.  IV. 6,  100  m  horizontal  resolution. 

Fig . 

IV.  8: 

Example  of  stream  line  wind  field  from  neutral 

LINCOM  (10  meters  height)  at  Porton  Down  ter¬ 
rain. 
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